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I. INTRODUCTION 

This document deaorlbss s computer package which produces contour plots 
of the three-dimensional electron distribution function measured by the 
GSFC eleotron spectrometer aboard ISEE-1 (Ogilvie et al. , 1978). The 
instrument configuration consists of three opposed pairs of analysers whose 
alignments are mutually orthogonal. Thus, at any Instant, the Instrument 
is looking along both senses of each of the three orthogonal axes, as shown 
in Figure la. The timing and en«*gy steps of the electron energy spectrum 
measurement is shown in Figure 1b. In one sweep through the 16 electron 
velocity steps of the instrument, the spacecraft rotates through one 
sector, or one-sixth of a complete revolution sampling phase space at 96 
(16 X 6) points in 0.5 seconds. The electron velocity sweep repeats during 
each sector and one complete 3 second spin of the spacecraft observes at 
576 (6 X 96) discrete points on the distribution function, f(^(t)), which, 
as a surface, resides in a 4-apace whose coordinates are (V^(t), V^(t) , 
V^(t), f(Vjj(t), V^(t), V^(t)). The implicit relationship on time given by 
^(.t) is due to the magnitude and direction of the sampled velocity vectors 
being determined by the energy sweep and the spin rate, respectively. 

A mod si- independent graphical representation of the measured f(V(t)) 
is illustrated in Figure 2 which shows 2-dimensional cuts of the full 
f(^(t)) surface for one spin. The two dimensional trace is defined by 

f(^ • ftj^(t)) vs ? ♦ <!^(t) > 0, and 

f(7 • fijtt)) vs ^ • nj(t) < 0 

f'tr opposed sensors which by hardware coastruction have anti-colinear 
fields of view, * '"j * opposed sensors have a common 

symbol. Thus each panel ir Figure 2 displays data from three pairs of 
opposed sensors rotated into a common plane to facilitate comparison. The 
six panels show data from successive 0.5 tecond energy sweeps and 60” 
rotations of the spacecraft. Although each "slice" of the f(V(t)) surface 
shown by these 2-dimensional traces allow model-independent determinations 
of aniso*‘ropies in the distribution function, a more complete picture of 
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the plaana microstate could be obtained if the four dimensional f(V(t)) 
surface itself could be represented. 

Except in high gradient regimes, the electron plaana is gyrotropic in 
the proper frame, l.e. , the frame that is at rest relative to the bulk 
flow. This results in a cylindrical synunetry about the magnetic field 
direction which reduces the 4-space to three-dimensions (V„(t), Wj(t), 
f(W«(t), Wj(t)), thus collapsing the distribution to reside in a 3-space 
which can be contoured. The electron speeds W„ and are the proper frame 
parallel and perpendicular components of the observed velocities. This 
transformation has been used to assemble the 2-dimensional cuts of the 
distribution function in Figure 2 as the contour plot of Figure 3a. The 
coordinates are and W„ and the contours indicate levels where the 
distribution function is a constant. The separation of the contours 
correspond, to variations in the natural logarithm of the distribution 
function equal to unity. The dotted lines show the "tracks" in Wj^, W„ 
space along *diich the observations actually lie, and the gaps between 
observations are filled, in this space, with a two dimensional spline 
procedure. In Figure 3b, a gray shaded version of the contour plot of 
f(W„(t), Wj^(t)) is shown. The level of gray shading is proportional to the 
logarithm of the distribution function and the white lines indicate 
discrete iso-contours. 

Either of these plots are a product of this computer package at the 
option of the user. The line contours provide a higher resolution display 
than the gray shaded contour u but require more computer time, approximately 
45 seconds CPU/IO time for each line plot and approximately 10 seconds for 
each gray plot. The remainder of this report will document the use of the 
program: Section II contains a description of the output plots; a brief 

user guide describing the user requirements and options are given in 
Section III; and an outline of the method of the synthesis analysis in 
Section IV. 
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II. Dtaerlptlon of tht Output Plots 
A. Lln< Plots of tn f lao-Contoura (Figure 4) . 

TWo rtprestntatlons of tho synthtalztd distributions are shown in Figure 4. 
In Figure 4a, the coordinates are the proper frame speed W and pitch angie 
e, where W a |8| * |v-u|, v • electron velocity, u * electron bulk velocity 
averaged over one spacecraft rotation (3 seconds) and 6 is the angle 
between W and B, the symmetry axis. In energy mode 1 or 2, B is along the 
averaged heat flux vector, and its sense is in the earthward hemisphere. 

A 

In mode 3, B is along the '.oCuetic field direction, obtained from the 
onboard magnetometer. In figure 4b, the coordinates are the proper speed 
components, and W„, and the plot is required to be symmetric about 

s 0. 

The solid lines are the iso-contours of tn f, successive contours 
corresponding to variations in tn f equal to unity. The dashed lines are 
the 36 "data tracks" in velocity space along which measurements are made 
during one spacecraft revolution (6 detectors x 6 angular sectors). In 
accordance with the gyrotropy assumption, all 36 data tracks have been 
rotated about the field direction into the same (W, e) plane and, 
therefore, some rotated tracks will appear to intersect in this coordinate 
system as shown in Figure 4. Since the electron plasna is sampled at 
discrete points in phase space, there will be regions in phase space in 
which there are no measured points which are filled in by interpolation in 
order to do the contouring. The Interpolation is done in two steps in the 
(W, e) representation; first, an Interpolation is done over the measured 
energies along each data track (l.e., each detector and sector) followed by 
an interpolation over pitch angle. The resulting surface in f as a 
function of (W, e) is then snoothed before contouring. The representation 
in (Wj, W„) coordinates in Figure 4b is obtained by a point-for-point 
transformation of the »oothed distribution in (W, e) coordinates in Figure 
4a. 
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B. Gray-^hid<d Plots (Figurt 5) . 

Th« plots in Figura 5a and 5b ara tha gray-shadad version of the line 
plots in Figura 4a and 4b. raspactivaly. Lina contouring provides batter 
pixel resolution, but gray shading requires less execution time. (The gray 
shading package used was prepared and documented by M. L. Kaiser and P. G. 
Harper.) The continuous gray scale is obtained by digitizing tn f to a 
scale of integers IZ given by 

IZ s (tnf - (tn f) max tn ZG) * NP/tn ZG 

where HP a 16 on 12" Versatec plots, and HP a 64 on 22" Versatec plots. 

The parameter NP, the dynamic range tn ZG, and the maximum, (tn f) max, of 
the tn f scale may be adjusted as discussed in Section III. The set of 
integers IZ become the element values of an input array which converts IZ 
to proportional shades of gray; IZ a NP is blackest, IZ a 0 is whitest. If 

tn f exceeds the allowed range, then IZ is set equal to the appropriate 

bound, IZ a NP or IZ a 0. For a given contrast or difference in gray 
levels, A(IZ), the corresponding difference in tn f is A(tn f) = 

(tn ZG/NP) 4 (IZ). 

The white contours are levels of constant tn f which are superimposed 
on the continuous gray scale. The fixed values of gray level at which 
contours are placed are the set of integers 

IZ (fixed) a IZ/(NP/8), IZ a 0 to NP. 

The difference in gray level between successive contours is A(IZ) ^ NP/8. 
The corresponding change in tn f is A(tn f) a tn ZG/8. For example, if the 
range of f(v) on the gray scale is set at ZG a 8.886 x 10^. then between 
fixed contours 4(in f) a (tn 8.886 x 10®)/8 a 16/8 a 2, or f changes by the 
factor e"^ between contours. 

The black bands running more or less vertically in the W -e grid of 
Figure 5a outline the data void in the parallel and anti>parallel 
directions. These are the locations in (W -9) space of measurements of the 
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two dotoetora orlentod eloaost to 6 ■ 0* and • ■ 180* pitch angle. The 
region between the outer banda and the edge of the grid, therefore, contain 
Interpolated data only. Aa will be dlacuaaed in Section IV, the 
Interpolation la performed by noting that f(W„,Wj^) a f(V„, -W^) , i.e. , that 
f ia mirror aymmetrio about the B direction. Thua the interpolation near 
0* and 180* ia not an extrapolation; rather it ia a trend ayntheaia. 

In Figure 5c, the parallel and perpendicular cuta of the anoothed and 
Interpolated to f are plotted aa indicated on the plot. Note that theae 
cuta may not be determined by direct meaaurementa in the aelected dirctlon, 
rather, they are the trend auggeated by the cloaeat meaaurementa and the 
"aperture ayntheaia" approach adopted here. 

Parallel: in f(W,e * 0*) ia taken from the W -e grid and plotted aa in 

f(W„) va on the right aide of the plot. On the JLeft-hand aide, to f 
(W,e * 180* ia plot‘;ed aa to f (-W„) va -W„. 

Perpendicular: in f (W,8 s w/2) ia alao taken from the W -6 grid and 

plotted aa to f (W^) va W^. The plot ia aymmetrlc about * 0 due to the 
required aymmetry about B. 

One-count level: The in f corre^onding to the one-count level of the 

detector ia plotted on the aame acale aa the parallel and perpendicular 
diatrlbution. It la not a anooth curve aince the integration time ia not 
independent of the aampled apeed. 

Vertical Axia: The acale of in f ia from (in f) max to (in f) max 

-in Z; (in f) max can be the aame value aa uaed in the contour plot a or it 
can be a aeparate fixed value aa dlacuaaed in Section III. The range of 
to f, i.e., to Z, ia a required input and la dlacuaaed in Section III. 

In Figure 5d, the coordlnatea axia are the aame aa thoae in Figure 5b. 
The bottom half plane contain a the velocity apace location a of the meaaured 
data polnta during one apacecraft revolution. Theae are the pointa uaed 
in determining the data tracka ahown in Figure tt. The upper half plane 
ahowa contour a of the ayntheaized diatrlbution function, but without any 
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anoothlng. This plot should bo comp#, sd with the anoothsjd plot in Figuro 
5b to gsin an improsslon of th# levol of fluctuations in tho synthosizod 
distribution ind what tho offoct of anoothing has boon. Tho anoothing usod 
aftor tho intorpclation is to proeludo tho introduction of angular 
signatures sharper than tho angular homogmoity of tho analysis system at 
any given «iorgy, viz., otr (4ir/6 detectors x 6 sectors) x 57 * 3*^ 
degrees. 

In Figuro 5o, tho natural logarithm of the reduced distribution 
function, in F(W„) (vortical axis) is plotted againt'it W„ (horizontal axis). 
The reduced distribution is given by 

m 

F(W„) s 2v / f(W„,Wj.) Wj d W^. 

0 

Tho program computes F(W„1 by Integrating the anoothed, un-norraalized 
values of W, f(W„,W.) from Figure 4b along columns, i.e. , over at given 
W*. The in F(W„) is plotted on a scale from in (5 x 10 ) to in 

(5 X 10**^^) , allowing comparisons to be made between various distributions. 
The horizontal scale 1s identical to that of Figure 4c. 

The Landau damping decranent of plaana oscillations that may be 
present, in (-X) is computed and plotted in Figure 5e as a function of W„. 
The scale of in (— x) is from tn '(10^) to in (10 ), Details of the 
calculation of F(W„) and x are given in the Appendix. 

Points for which x > 0 are ignored at present. This of course is the 
interesting case in which a positive slope of F(W„) results in the growth 
of plaana oscillations. In a future version of the progran, these points 
could be flagged and included on the same plot as the damping decrem«it. 

The two printed lines at the top of the plot in Figure 5 are selected 
electron fluid moment parameters averaged over one spin (3 seconds). 

First line: N * electron density (cm '') 

U s electron bulk speed (km/s) and its 
GSE X,Y, Z components 
T s temperature ("K) 

A s anistropy, T„/T^. 
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Second line: H s heat flux (erga/cm^a) and ita 

GSE • coordinatea 

S c algn of ® *^E 

IT 

Wien S « -1, fl polnta toward Earth. 

S 3 <^1, 0 polnta toward the Sun. 

The date when the plot waa made la at the far right at the top. At the 
bottom of the plot, the year, decimal day, and UT of the meaaurement la 
printed. 


Ill* Ualng the Contouring Program 


A. Required Input a 

The required inputa are the aame for the gray-ahadlng and line contour 
veraiona and are given below: 

JDECDT 3 decimal day of year 

JHR 3 hour a Day end time (UT) of the atart 

JMIN 3 mlnutea of the data interval 

JSEC 3 aeconda to be proceaaed. 

MAPS 3 total nwBber of apectra (one each spin) to be 
plotted. 

ISKIP 3 a parameter determining whether conaecutive apectra 
are to be plotted. 

For example: ISKIP 3 1, every apectra la plotted; 

ISKIP 3 2, every other one, etc. 

The time in aeconda apanned by the data Interval ia MAPS»ISKIP * 9 

or MAPS»ISKIP*18 depending on the data format: In format 1, the apectra 

are 9 aeconda apart; In format 2, the apectra are 18 aeconda apart. 
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B. Program Options for tht Gray-Shading V»r»l_on . 

Various parametars which control tha availabla options ara antarad 
through NAMBUIST/INPUT/for tha gray varalon; thara ara no options for tha 
lina oontov. xng var sion. 


1. Plot parametars 


Paramatar 


Function 


Default 


JPRNT a 0 Various moments computed over in- 

dividual sectors and ^In-averagsd ara 
Included in printout. 

s 1 Extended printout lists input data and 

tabulation of_contourad data. 


IZSET a 0 

a 1 


ZMXX 


Tha f*,fj plots use default value of 
scale maximum. 

ZMXX will override default value 
(Note; If IZSET a 1, value of ZMXX 
must be entered). 

Desired scale maxlravro of f„.f^ plot. 


0 


1 x 10 


-24 


IZABS a 0 The f«.f^ plots are normalized to 1. 

8 1 A fixed scale maximum of fo*" ■H 

plots will be used. The value used 
depends on IZSET and ZMXX entered above. 


IFLOW a 0 


Plots assume proper frame correction is 
negligible. 


a 1 


Z 


Plots are in proper frame coordinates. 
(•Note; When IMODE a 3, default a 0). 

Dynamic range of f on f„. f^ plots 
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20 Dynnlo rang* cf f oovtrlng rang* of 

gray acal*. 


1 X 10^ 


paramatara ZMXX« Z and ZG ar* R*X variablaa; the rcat ar* l*t. 

The ‘ollotelng axamplaa illuatrat* aom* poaalbl* plot optiona: 

a) Change nothing: If only default values ar* used, then ZG s 1 x lO^t 

2 ■ 1 X 10^ and (f) nax on f,,fj^ plot will be 5 x 10"^^ for modes 1 

and 2, and 3 X 10“^^ for mode 3. 

b) To specify rang* of f on contour plot: change ZG. 

c) To ^eclfy rang* of f on f„, f^ plot: change Z. 

o) To specify f max on f^.f^ plot: set IZSET « 1 and specify ZMXX. 

*) To aorraalixe f„,f^ plot to 1: set IZABS * 0. 

2. Additional parmneters 

The following parameters controlling output have been added to 
NAMELIST/INPUT/: 

Parameter Function 

KPRNT d 0 The values of tn f on the (W,e) grid may be 

printed out at various steps In the syntneals 
analysis. (See the current listing for details.) 

Default: KPRNT s 0 

ITAP ^ 0 The reduced distribution function as a function 

of parallel proper velocity is written out on tap*. 


Default: ITAP * 0 
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3* Synmttry axis 


In tha currant varaion of tha gray^ahadini package, the aymmatry axla 
with unit vactor componanta B(I), 1 * 1.3 <CSE Cartaaian coordlnataa) may 
ba takan to ba along althar 1) tha alactron haat flux vactor direction with 
Ita Sanaa earthward, or 2) along another apacifiad direction parallel to 
tha vactor with componanta BF(I). 1 « 1.3* The parameter controlling tha 
uaara option la IBEXT and la entered through NAMELIST/IHPUT/ . The option a 

are a a follows: 


Paramatar 


Function 


IMODE * 1 or 2 
IBEXT s 0 


Tha aymmatry vactor B la along tha alactron 
heat flux vactor but in the earthward sense. 


IHODE « 3 
IBEXT a 0 


Tha aymmatry vactor B la along tha normalixad 
BF direction, where BF la tha 5 min average of 
tha magnetic field obtained from tha ISEE-1 
magnetic field experiment and is read from the 
input tape to this program along with tha 
alactron plasna data. 


IMODE « 1,2, or 3 Any symmetry vactor could ba apacifiad by 

IBEXT i 0 entering BF(I) I . 1,3. t.irough an additional 

input which overrides tha input magnetic 
field whan in IMODE . 3. The repulrad 
additional input is not now a part of tha 
program, but tha nacasaary branching for this 
option is in place. 


Default: IBEXT * 0 


tt. S.aoothing Parameters 


Tha rmsalnlng control parameters that may ba changed tnro'jgh 
NAMELIST/IHPUT/ are IS2T, NC. NH. WC. WH. These parameters control tha 
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liz* of tht anoothlng Inttrval in pitoh mflt and apood. It la not 
ooooasory to ontor than as inputs slnea thay hava dafault valuas. Ihasa 
parssatars and tha aa'. jthlng art daaorlbad in Sactlon ZV-E. 

B. Job Subwittal 


Tha followinj aaquanea of JCL statananta and data ara raquirad to run 
tha program on tha 360/ 91. Oaoks ara availabla on raquost. 

1. Cray-ahadlng varalon . 


00100 // EXEC LOADER. PE0I0N«499K.PAR«-'SI2E-495iC.EP*MA INIS. LET^ 

00110 //SYSLIB DD DSN*t12.U3MLK.L0DI1OD.DISP»SHR 
00120 // DO DSN-SYSl.f1002.F0RTLlB.DI£‘R*SHR 

00130 //SYSLIN DO 

00140 // DD 0SN-SYS2.L0ADLIB(U3RJFCRl).DISP*ShR 

00141 //OO.FT06OO1 DD SYS0UT«8.DCB«(RECFM»VBA.LRECL«137.BLKSi2E-l7S9> 

00142 //OO.FT08FOO1 DD DUMMY.SYS0UT“A.DCB*<RECFM»VBA.LRECL«137.BLKSI2E»1789) 
00130 //OO.FTlOFOOl DD DSN«ISEEA.SUMWARY.LABEL-( .SL. . IN) . 

00160 ft DCB*<RECFM-VS.BLKSIZE*14C56.LRECL«14052.DEN-4.0PTCD«B.BUFN0»1) . 
00170 / / V0L-SER-FP8804 . D I SP» ( SHR . KEEP > . UN I T« ( 6250 . . OEFfcR ) 

00180 //00.FT20F001 DD UNIT-800.LABEL«( l.NL. .OUT) . 

00190 / / DSN-REDUCE . DCB- < RECFM-F , BLKS 1 2E-4745 . LRECL-4745 . DENr*2 . OPTCD-Q ) . 
00200 ft D I SP» ( MOD . KEEP ) . V0L«SER»FP872 1 
00210 //DATA5 DO » 

00220 tc I NPUT I Z ABS- 1.Z-1E7.Z0-1E7.I ZSET- l.ZMXX-lE-24. JPRNT «0 . KPRNT-0 . JtEND 

00221 311 02 11 00 01 01 

00280 //VPLOT DD UN IT-800. DrD»DEN»2. LABEL* ( .NL) . 

00290 // V0L-SER-TD5054 
00300 // EXcC NTSO 


Linas 150>170: 
Lines 180-200: 
Lina 220: 

Lina 221: 

Lines 280-290: 


Input tape DD st<«tament. 

Optional output tape DD statauant (See Section III B.2) . 
NAMELZST/INPUT/ containing parameters controlling 
program option. 

Data Interval JDECDY. JHR, JMIN, JSEC. MAPS. ISKIP In 
the format 13, 5(1X,I2). 

Output plot tape DD statement. 


Approximate execution time: 10 sec CPU 


1 sec 10 
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2 . Line plot v«*sion. 


OOOlO // EXEC LINKG0»REGI0N-499K.PARM«'LET' 

00020 //SYSLIB DD DSN«SYS2.U0LF. VERSATEC. DISP«SHR 
00030 // DD DSN»SYS2.W0LFPL0T,DISP»SHR 

00035 // DD DSN-SYS2.IMSLS.DISP-SHR 

00036 // DD DSN-SYSl.FORTLIB.DISP-SHR 

00037 // DD DSN-SYS2.F0RTLIB.DISP-SHR 

00040 //OBJECT DD « 

00050 INCLUDE LO ADLIB (U3RJFL IN) 

00060 INCLUDE SYSLIB( lEVMAP) 

00070 ENTRY DFCN2R 

OOOSO //GO.FT08FOOI DD DUMMY.SYS0UT«A,DCB»(RECFM-»VBA,LRECL»137,BLKSIZE«1789) 
00090 / /CO . FT I OFOO 1 DD DSN» I SEEA . SJJMMAR Y » LABEL- ( . SL . . I N ) , 

00 1 00 / / DCB- ( RECFM- VS » BLKS I ZE- 1 4056 , LRECL- 1 4052 . DEN-4 , OPTCD-B » BUFNO- 1 ) . 

00 1 1 0 / / V0L-SER-FP7877 , D I SP» ( SHR , KEEP ) . UN I T» ( 6250 ♦ . DEFER > 

00120 //PLOTLOC DO SYSOUT-A 

00130 / / VECTR I DD DSN-&J* VECTR 1 . UN I T-SYSDA , D I SP- ( ♦ PASS > , SPACE- ( TRK . ( 1 , 1 ) ) 
00140 / / VECTR2 DD DSN»«<StVECTR2 , UN I T-SYSDA » D I SP- ( . PASS ) » SPACE- ( CYL . tl . 1 ) ) 
00150 //0ATA5 DD ♦ - . 

00160 &INPUT WX»1...5»l.,.5tFLD»2..2..l..l..ScEND 
00170 312 12 33 00 01 01 

00180 // EXEC PCM-IEVMAPP.COND»(EVEN),RECION-300K 
00190 //STEPLIB DD DSN-SYS2. WOLF. VERSATEC.DISP-SHR 
00200 //PLOTLOC DD SYSOUT-A 

00210 / /VECTR 1 DD D ISP- (OLD. DELETE) ♦ DSN-Sii VECTR 1 
00220 //VECTR2 DD DISP*(OLD. DELETE) . DSN-&5cVECTR2 
00230 //RASTTAPE DD UNIT-800. DCB-DEN-2, LABEL- (. NL) . 

00240 // VOL-SER-TD5060 
00250 // EXEC NTSO 


Lines 90-110: Input tape DD statesont. 

Line 170: Data interval JDECDY, JHR, JMIN, JSEC, MAP?, ISKIP in 

the format 13, 5( IX, 12). 

Lines 230, 240: Output plot tape DD statement. 


Approximate execution time.- 45 sec CPU, 1 sec 10 


ZV. Method 


This section outlines the method of synthesizing the discrete 
measurements of the distribution function In three dimensions into an array 
of continuous values which can be contoured. Only the essential steps of 
the analysis are outlined here; the details are found In the docmented 
computer code which Is available upon request. The parts of the computer 
code that are relevant to this outline can be Identified by commrnits in the 
code corresponding to the subheadings found here. 

A. Read Input data . 

1. Read the control parameters (explained in Section III). 

2. Read the header record containing physical characteristics of the 
detector which are used In computing the one-count phase density for the 
given mode of detector. 

3. Read the start time of the data Interval and the number of records 
to be processed. 

4. Read the data records in the data Incerval and process one record 
at a time. The essential data that are read are: 


IYR,IDECDY.IHR,MIN,SEC 
CC(J), J s 1, INSET 


CMTS< J.L.IS) 
FBLOK(J,L,IS) 


- Date and time of start of data interval, 

- Array of ^eeds corresponding to the 

tuned energy of spectra, in order of 
decreasing magnitude, 

- Electron counts, 

- Phase density of electrons. 


where J is the velocity index and L,IS are the detector and sector 
Indices, respectively. 


IMODE 

SPINP 


- Energy mode. 

- Spin period 


Various moments of the distribution function computed for each 
sector and averaged over spin are also read. 
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B. Computt thf atimuth In aoln plan* of tht lo ok dlrtctlon of t«eh 
dptfotor In tach apctor . 

PHIZ(L.IS) s (360/SPlNP)«CYCLT»(IS-1)^DPHI(U±90 for IS * 1, 6, 

wh*r« SPINP « upditPd spin p«riod 

CYCLT « atctor timt 

0PHI(L)±90 art the datfctor azimuths ralatlv# to tha saetor and since 
DPHId) * DPH(6), DPMK2) « DPHK5) . and DPHK3) « DPHK4), tha 
azimuths DPHI(L) t oq* ara 180* apart for tha datactors in aach of the 
pairs (1,6), (2,5), and (3.“)* For L * 1,2, or 3. the lag is -OO*; 
for L » 4,5,6, the lag is +90*. 

C. Choose the axis of aymmatrz » 

• 

Whan IMODE * 1 or 2, tha symmetry axis fi is along tha heat flux 
direction H. The sense of 6, however, is chosen always to be directed 
toward the earthward hemi^here. Thorefore when H points toward the earth 
and away from the sun, fi « fi/H; when fi points away from the earth, fi t 
-H/H. Wh« IMODE « 3, the symmetry axis must be provided independently of 
the electron data. (See Section III.B.2). 

D. Transform velocity space from the S/C frame to the proper frame of 
the electrons with 6 as the axis o f symmetry. 

1. Form the array V(I),I « 1.NV. of the measured speeds from the 
input array CC. The V array is ordered according to decreasing magnitude. 

2. Order the measured speed array, V, and sample time array, TIM, 
according to increasing velocity magnitude and name then VS(I) and Y(I), 
respectively. Compute the Interpolating spline over Y(tlme) as a function 
of VS (velocity magnitude). See Figure 6a. 


3. aear the plotting arrays. 


16 


4. For givon dotoctor ^nd sector indieest L snd IS, eliminate 
array elements having zero or negative (l.e. . "fill") electron counts and 
form the new arrays 

TP(I) « time 

YP(I) « log e of the distribution function (in f) 

VP(I) s measured velocity magnitude 

I s 1.NV 

These arrays are ordered according to decreasing velocity 

magnitude. 

5. Order the arrays TP.YP, and VP according to increasing speed, 
forming the new arrays 

T(I) X time 
B2(I) s tn f 

Bid) « measured veleocity magnitude 
1 s 1 ,NV 

Compute the Interpolating spline over 62(in f) vs. B1 (speed) 
for given detector and sector corresponding to Indices L and IS, 
respectively. See Figure 6b. 

6. The transformation and interpolation 

In order to contour in f measuronents, the data are 
transformed to a regularly spaced grid in the proper frame coordinates W,e 
where W is the proper frame speed and o is the pitch angle relative to the 
symmetry axis S. Ihia is done by choosing values of W and, by iteration, 
finding values of (V(t). satisfy 

(w')^ X 

X y z 
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which, using V detected « V detected bi*/ be written in componenta 

aa 

Wjj « -V(t) ain6sg(t) coadggCt) 

Wy s -V(t) ainesg(t) aindjgCt) -Uy 

Wj s -V(t) coaiggCt) -Uj 

(Notice that a, a are the "look" angles rather than the acceptance angles 

of the collected flux ■ “Kdateeted*^ Since V(t) and («(t), a(t)) 5 g 

depend on one another because of the timing sequence through the energy 
channels and the spin rate of the S/C, values of V and («,a)££ used in the 
iterative process, must be consistent with the functional relation shown in 
Figure 6a for V(t) and with _ 

Sgg(t) s THETA(L) s const for given L 
agg(t) * PHIZ(L,IS) ♦ (360/SPINP)»t. 

Using the computed spline over V(t) in Figure 6a, the value of t s t* 
and hence V s V that yields the desired W* is found, and the pitch angle 
e* is then given oy 


« cos 


-1 


''x ®x " ^ ®y ^ ''z ®y 


Ambiguity in cos"^ is irrelevant because of gyrotropy and, therefore, the 
principle range is 0 to t. The correspondilng value of in f* at V s V is 
then found using the spline over in f vs V given in Figure 6b. 


The reajlt of transforming the velocity coordinates of the measurement 
from the single detector and sector given in Figure 6b is shown in Figure 
6c. This data "track" has data points at every horizontal grid line 
(values of VO which in general, don't fall at the intersections of the 
vertical grid lines (values of 6p). After all 36 data tracks have been 
transformed, intersections of all the tracks with a given horizontal grid 


line arc uaed to compute an interpolating apline to find tn f at each of 
the 121 grid points on the glv«i horizontal line. Figure 6 d illustrates 
this horizontal interpolation. The mtire grid is filled by repeating the 
procedure for all 60 horizontal grid lines. The interpolation is performed 
in subroutine INTRP 3 . This completes the synthesis of the measurements in 
the (W,e) space. 

A cautionary comment is in order regarding the horizontal 
interpolation. Note in Figure 6d that the intersections of the 36 data 
tracks with any given horizontal line are unevenly spaced, being widely 
separated in some places and crowded or overlapping in others. Since the 
interpolating spline fitting in f vs 6^ will pass through every point on 
the line, those points that tend to overlap in 6 may cause the spline to 
ring. This can lead to spurious interpolated values depending on how large 
the gap la between the overlapping points and the nearest neighboring 
points along the line.* The approach to this problem is to average together 
the points in a given row into 6 degree pitch angle bins before doing the 
interpolation. This Irreducible bln size is fundamental and determined by 
the field of view of the sensors which is 8.5 x 11 degrees in solid angle. 
The synthesis therefore reflects solid angle averages which this finite 
solid angle implies. 

The data gap which is usually largest is at the edges of the W -e grid 
near e = 0* and e s 180* , resulting from the fact that the magnetic field 
is usually close to being in the ecliptic plane when in the solar wind 
while no detector views the sun directly. To fill this gap with 
interpolated data and to minimize ringing effects, the spline begins and 
Olds outside the space of Interest, i.e. , outside the space e = o” to 180^. 
This is done by reflecting the half-space between 9 s 90 * to O" about 9 s 
0° and reflecting the half-space between 9 s 90 * to 180° about 9 s 180”. 
This is consistent with the assumption of cylindrical symmetry. Therefore, 
tr.ese gaps are filled with Interpolated data and are not an extrapolation. 
Note that in the absence of any tendency for iso-contours of tnf to slope, 
then the de facto interpolation across the data gaps at O” and I 80 ” will be 
contours which are straight lines in (W,9) space and circles in (Wj^.W„) 
apace. Therefore if the synthesis suggests non-circular deformation at 0“ 
and 180”, it must reflect contours with 3f/39 iO. 
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Smoothing the Inttrpolattd data . 

Before contouring the in f date on the filled W -e grid. . two variable 
anoothlng function S(W,«) la uaed. The characteriatica of S(W,e) , which 
slldea over the grid, are aketched in Figure 7. Since the channel 
eeparatlon varlea with apeed. the width along W ia variable, being wider at 
higher apeeda. The parametera ISZY.WC.WH.NC. and MH are adjuatable through 
NAMELIST (aee Section IV) in the gray-ahaded veraion. In the line plot 
veralon. theae parametera are fixed and equal to the default vaiuea of the 
gray— ahaded veralon. The default vaiuea are: 

ISZY a 1. (a R*4 variable) 

NC « 10 unita of pitch angle (1 unit a 1.5 degreea) a 15 degreea 

NH a NC 

NCal.jWHaO.^ ~ 

Thu, th. core of the o.f.uU I, 30" old. .„d 1, n.t lo «putude; 

the wing a extend 15* on either aide with reduced amplitude. 

Thia choice of anoothing function waa baaed on the following 
conaiderationa: 

a) Coaraeneaa in the data ia introduced by the two-way interpolation 
aa dlacuaaed above. 

b) The aize of the anoothing window in the 0 dimenalon ahould not be 

finer than t.he degree of continuity that ia Implied by aamplir.g f(V) in 36 
diacrete directiona over 4, ateradiana at a given energy during one 
apacecraft rotation. An eatimate for that average apacing in angle ia v* 
(ttt/36) • 57 degreea s 3U degreea. 

c) The computation muat be economical in computer time. Thia ia why 
the anoothing function waa conatructed uaing flat aegmenta. 
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F. Transformation from (W, e) to (W,, %) eoordlnatts . 

Tha anoothad data on the (W, e) grid are transformed to the (W„,W^) 
grid with no additional anoothlng. A two-dimensional Interpolation Is used 
as followa: 

G’’ X in f (W^.Wj) r 

(W'-W^) [(e* -e^) in f («j ^ ^ 4 . 1 ^ ^®1*'^1 ♦ 1^^ 

(Wi ^ ^-W) C(e' - e^) in f (s^ - «'> ^ (e^.w^)] 

VIhere W„,W, Is a point at an Intersection of grid lines, W',e' Is the 

2 2 

corresponding point on the (W, •) grid, W * , and 6* s 

cos“\w„/W'). The constant G s (Wj^ ^ ^ - Wj^) ^ ^ -d). 

The point (W, e‘) 1s not at a grid Intersection, In general, but falls 
within the box whose corners are (Wj^,e^), (Wj^ ^ ^ ^) , and 

(Wi ^ i« ®1 + 1^* Interpolated value contains contributions from the 

value of in f at all four corners weighted by the Inverse of the distance 
of the point (H*,e*) from each corner. 

The full polar plot In Figure 4 and Figure 5 Is obtained by reflecting 
the plot about s 0, l.e., 

in f (W„,-W^) s in f (W„. Wj). 
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FIGURE CAPTIOUS 

ES“SE_1 (.) V.locltif «P«. of tho (b) Ttalo, 

diair. for tlio ,l«.tron «>«■,, mnmirmmt. 

HSIBLl th. ..loolt, dl«rlbotlo» functloo bf tb. 

ISEE VES burloj Ob. of th. 

. ap«c#er«ft. lh« data in aach sector is acquired In 0.5 
*«.bd.. th. thr.. p.ir. Of oppo,«l «„»r. b.ln, »,pi« 
almultMMuair in .Mh ««-,y ch.nn.l. Ih... 576 dlKr.t. 
■««.r«„t. T. dlarlbot« throuphout ph.a a.o. du. to 
th. tlmlb, of th. «a,y «„.. ,0 th. v.c«r.ft (o.d throuph 
one complete rotation. 

liSHSLl (■) Line-contour plot of the synthesis of the discrete 

«-ourm„ents of the distribution function shown in Figure 2. 

velocities, and W,, are the speeds perpendicular and 
P-rallel to the magnetic field in the frame moving with the 
« ec ron bulk How. Cylindrical symmetry about the parallel 
direction Is assumed, (b) Cray-shaded version of contour 
^ot Showing gradual changes in the distribution function, 
e white lines are contours at fixed values of f(tJ). 

IISliSLl Output of the line contour plot package. See text. 

miHLl Output Of the gray-shaded plot package. See text for 
detailed explanation. 

£J£!15S_S (•) Tlbin, Mqumc. of th. v.locltl.s it which th. 

dlitrlbutlob funotloo Id umplri ,t ..oh d.t«tor In ..oh 
d«:tor. (b) „lu., of lo, . of th. dldtrlbutlon 

function, to f. ,t . p.rtloul.r d.t«tor (L . 1 ) .„d ,«:tor 
(IS . 1). By flttln, apiin.d to th. dati In (.) .nd (b) 
point, ir. fUl.d in by lnt«-pol.tlon thl, phii. ip.o. 
trick lod r.pr.d«,t« 1 „ propir friw, ooordln.t. ip.M i„d 
Pitch in,l., I, In ( 0 , Tb. looitlon of lU 36 d.t. trick. 
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during a oomplata rotation art ahown In (d). Tht pointa 
along aaeh track art ehoaan to Intaraaot horizontal grid 
llnaa, but, In ganaral do not Intaraaot at vartloal grid 
llnaa* Tha valuaa of in f and pltoh angla at tha horizontal 
grid intaraaotlona (ahown for ona horizontal grid llna) ara 
uaad to Intarpolata at tha intarvanlng vartloal grid 
Intaraaotlona. Tha grid aiza la 60 (apaad) by 120 (pitch 
angla) . 

FIGURE 7 The two«dimenalonal anoothlng function that alldaa over the 
V,a grid, anoothing tha valuaa of in f filled In by tha two 
way Interpolation (aaa text). 


APPENDIX: Calculation of tht r<ductd dlitrlbutlon fun ction and tht 

Landau damplna dacraaant. 
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. Tha raducad diatrlbutlon function lf(W,) that la plotted In FUura 5e la 
defined aa 

* 

F(W„) « 2» / dWj^ f(W„, Vj} 

0 

where and W, are the parallel and perpendicular proper frame velocitea, 
reapectlvely, and f(W,, W^) lathe anoothed, ayntheaized gyrotroplc 
diatrlbutlon functlona in the W„. velocity apace aa ahown In Figure 5b. 
The preclae llnita of Integration can extend only over the range of aampled 
data, to Hivever, the contribution to the Integral outalde the 

aampled"range of velocltiea can be approximated with help of the following 
change of variable: 


X « 


Y * 


- '^min 




min 


A2 


<«th ♦ '"a - «-ln^‘ 


f(W„, W^), 


"th 


vihere s la eatlmated thermal apeed. Then, f dW^ f(W„, W_^) » 
cn u 

^cr»ln 


/ Y dX. 


Extending the lower Integration limit to zero, » 0, can be done by 
noting that the Integrand vanlahea, l.e. , Y » f(W„, » 0 at * 0. 


The upper integration limit can be extended to « • by noting that 
* 1 aa W and aasuming that f(W„, W^) falla off faater than . 

S8X 

Therefore 
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F(W,) m 2« /' dx Y. A3 

0 

Th« varlabl«a X and Y ara glvan by Equation A2 and tha Intagral ia 
oaleulatad by an n-point trapaxoidal fomula. lha aat of (Xj. Y^) for 1 a 
2. n-1 ara takan from polnta along a vartlcal grid llna at givan v, in tha 
rapraaantation of Figura kb or Flgura 5b. Tha and points ara X^ * 

0, Y^ a 0 and Xjj a 1, Y^ a 0, eorraapondlng to a 0 and a •, 
raapaotlvaly. 

Tha Landau damping daeramant. W^, ia eomputad aa followa (Krall and 
Thivalpiaoa, paga 520): 


-1 



— 

^ aw. 



Ak 


*'*’*‘* “p* • alaetron plaaa fraquancy, k a wsva numbar, and a raal part 
of tha alaetron plaana wava fraquancy. 

If wa aat /k a and a a 2t 10** (n^)*^^, wnara n^ a electron 
density (ca~^)t than. 


f 

2 



a inF 


? k 2 

a - 10 n^ w,* r a inF/aw, 

Equation A5 ia used to eonput'a and plot in Figura 5a the plaana wave 
dmnplng daermnent as a function of W,. The plot of for all W, does 
not imply that plaana waves ara nacaaaarlly preswt at all wave velocities, 
but rather, that if waves ara present, then the damping effect of this 
electron distribution would be given by Equation A5. Polnta ara excluded 
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from tho plot which oorrcapond to wov# (rowth« l.t., when I in F/iW„ > 0. 
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